Scattering for the Gross- Pitaevskii equation 



Stephen Gustafson, Kenji Nakanishi, Tai-Peng Tsai 
February 2, 2008 

Abstract 

We investigate the asymptotic behavior at time infinity of solutions close to a non- 
zero constant equilibrium for the Gross-Pitaevskii (or Ginzburg-Landau-Schrodinger) 
equation. We prove that, in dimensions larger than 3, small perturbations can be 
approximated at time infinity by the linearized evolution, and the wave operators are 
homeomorphic around in certain Sobolev spaces. 

1 Introduction 

The cubic Schrodinger equation for ip{t, x) : M x M'^ ^ C, 

idtip = -Aip + \ip\'^ip, (1.1) 

has numerous physical apphcations. Noting that the single mode uniform oscillation e~** is 
trivially a solution, and writing (p = e~'^^ij), we arrive at 

= -A^ + (IVP - l)^', V(0,-)=V'o, (1.2) 

the Ginzhurg- Landau- Schrodinger or Gross-Pitaevskii equation, which models the dynamics 
of Bose-Einstein condensates, superfluids, and superconductors (see (3 for a review). In 
these physical settings, (p = e~** (or ^ = 1) corresponds to a stationary, constant-density 
condensate. Typically, one is interested in local perturbations of such an equilibrium, in 
which case the natural boundary condition as \x\ oo should not be the condition ip ^ 
usually imposed when studying (ll.lj) . but rather 

IV'I — > 1 as |x| — > oo. (1.3) 

Equation (|1.2j) with boundary condition (|1.3|) possesses, in dimension d = 2, well- 
known static, spatially localized, topologically non-trivial vortex solutions of the form 
^l^{x) = e*"''"''^^^^/(|x|), n G Z. There are also traveling-wave solutions, ip{t,x) = v{xi — 
ct,X2, ■ ■ ■ ,Xd), with V resembling a pair of vortices {d = 2), or a ring of vorticity (d > 3): 
[SI m IS For various Ginzburg-Landau Schrodinger and wave-type equations (such 

as p.2j) ). the dynamics of vortex points {d = 2), filaments {d = 3), etc., has been exten- 
sively studied in certain asymptotic limits: see, e.g., ^2 El HOI El UHl E] and references 
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therein. (In addition, there is a large hterature on vortex dynamics in parabohc versions 
of (HSl).) 

Very httle, however, seems to be known about the long-time behaviour of solutions 
of H1.2() - (|1.3() . and, in particular, about the stability of the various interesting vortex-type 
solutions mentioned above (however see, eg., I22j for numerical and formal asymptotic 
studies of vortex traveling wave stability, and ^21 for vortex stability /instability in the 
gauge-invariant wave-equation version of (|1.2|1 ). 

As a first step in this direction, we consider here a much simpler problem: asymptotic 
stability of the vacuum solution 

tpit, x) = e*^» = constant E S^. (1.4) 

By the gauge change — > e~*^°^, we may assume Oq = 0. Further writing tp = 1 + u, we 
have 

idfU + Au — 2 Re n = F(u), 

(15) 

F{u) := {u + 2u + \u\'^)u. 

The main goal of this paper is to investigate the asymptotic stability of the vacuum in 
terms of the scattering theory around the constant solution (|1.4j) of (|1.2j) , or the solution 
of TT^i . 

The first step is to introduce the following change of variable in order to represent the 
free (linear) evolution as a unitary group: 

v:=V-^u ■.= U-^Reu + ilmu, C/ := y^A(2^^A)^. (1.6) 

Then v satisfies the equation 

idtv - ^/-A{2-A)v = -iV-^iF{Vv). (1.7) 

The singularity of at the Fourier origin seems a natural effect of the interaction with 
the nonzero constant background in ()1.2|) . Due to this singularity, together with a quadratic 
nonlinearity which is not gauge invariant, the scattering problem for (|1.5|) turns out to be 
much more difficult than the original NLS (|1.1|) with spatially decaying data (i.e. (p G L^). 
Thus we have obtained so far the scattering only for higher dimensions d > 4. Global 
wellposedness of p.5|) in the energy space u G H^{M.'^) was proved in 3^ for d = 2,3. It is 
based on a priori bounds of the energy and so does not work for d > 4, where the cubic non- 
linearity is either energy-critical or super critical. Our theorem gives global wellposedness 
for d > 4 and small data by-product. 

Theorem 1.1 (Stability and scattering) Suppose that d > 4 and \a\ < — ^. // 

W^V-^uiO) is sufficiently small in ^'^/^-^(M^), then VV'^uit) remains small in H'^/^-^{R'^) 
for all t G M. Moreover, there exist f± G U^'^ H'^^'^'^ (W^) such that 

\\U^(V-Mt) - e-'^*i;±)||^./2-i(i,.) - (1.8) 

as t ^ iboo, where H = ^J—A{2 — A), and the wave operators v± ^ V~^u{Qi) are local 
homeomorphisms around in U~'^ H'^/'^~^ (W^) . 
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The corresponding result for the cubic NLS Hl.lj) holds |5] for any d> 2 and \a\ < d/2 — 1. 
Moreover, the wave operators are defined on the whole energy space for d = 2,3, 4, and 
are bijective for d = 3, 4 For d = 1, modified wave operators have been constructed 

in [23] for small data in weighted Sobolev spaces. 

The same result for NLS with the nonlinearity F(u) is easily obtained for d > 4 and 
\a\ < d/2 — 2 hj the argument in [31 . The lower dimensional case is more difficult due to 
the quadratic terms and slow decay of the free evolution. For d = 3, the wave operators 
have been constructed together with their inverses for the nonlinearities \u\u, and it} for 
small data in weighted Sobolev spaces [Hl^^- ^'^'^ d = 2, the modified wave operators for 
\u\u and the ordinary ones for u'^ and v? have been constructed for small data in weighted 
Sobolev spaces |H1E]- It has been proved that \u\'^ in d = 2 does not admit the ordinary 
wave operators ^5] . 

The estimates in weighted spaces are quite sensitive to the linear part and so do not 
easily extend to our problem. Here the difficulty in our equation is that the Galilean 
invariance is also destroyed by the boundary condition. For the Maxwell-Higgs equation, 
a similar hyperbolic problem, the Lorentz invariance is preserved and the scattering result 
has been obtained in j27j . 

To prove the above theorem, we first derive decay and Strichartz estimates for the 
linear evolution e~*^*. These are slightly better for d > 3 around the Fourier origin than 
the free Schrodinger equation, because of curvature properties of the characteristic surfaces. 
This is the reason why we have a wider range of a than the NLS with the same nonlinearity. 
Actually, this gain is strong enough for d > 5 to cancel the U^^ singularity in (|1.7j) . but not 
enough for d < 4. Instead, we use a seemingly miraculous change of the unknown function 
- a kind of "low frequency normal form" transformation - which completely removes the 
singularities at the Fourier origin from the nonlinearity. However, it still leaves regular 
quadratic nonlinearities, whose decay is out of our control for d < 3. 

In the next section 2, we investigate the linearized operator, deriving the decay and 
Strichartz estimates for the evolution. In Section 3, we perform the change of unknown 
function, and then in Section 4, prove the main theorem. 

2 Linearized operator 

In this section, we derive the decay and Strichartz estimate for the linear equation 



iu + Au — 2 Re li = /. 



(2.1) 



As in the introduction, we define the following C-linear Fourier multipliers 



U := V-A(2- A)-i, H := V-A(2 - A) 



(2.2) 



and R-linear operator 



Vu := U Re u + i Im u. 



(2.3) 
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By using the identities H = —AU ^ = {2 — A)U, we obtain 

iu + Au-2Reu = iV{dt + iH)V-^u. (2.4) 

Thus the hnear evolution by (|2.1jl is reduced to the unitary group e~'^^^, for which we have 
the following decay and Strichartz estimates. Here := denotes the homogeneous 

Besov spaces (cf. ^|), and we denote spacetime norms by := ^ II /(^) Ollx*^*^ 

etc. 

Theorem 2.1 Let d > 1. (i) For any 2 < q < oo, we have 



(2.5) 



where q' = q/{q — 1) is the dual exponent and a = 1/2 — 1/ q. 

(a) For j = 1, 2, let 2 < Pj, qj < oo, 2/pj + d/qj = d/2 and Sj 
{Pj^Qj) / (2,oo). Then we have 



^{1/2 - l/qj), but 



t 



11 

-i{t-s)H 



f{s)ds 



— oo 



(2.6) 



LPiBO^ 



Lp'2B0, 

"2 



When d > 2, these estimates are slightly better than those for the free Schrodinger equation 
because of the factor U, which is small around the Fourier origin. This gain is due to the 
fact that the curvature of the characteristic surface |CI v^ + lCP is larger than that of |^p 
in the {d — 1) angular directions and smaller in the radial directions. The above theorem 
follows from the standard argument for the stationary phase. For completeness, we give a 
proof in the following subsection. 



2.1 Stationary phase estimate 

Theorem 2.2 Let ip{r) G C°°(0, oo) satisfy the following. 

1. (p'{r), ip"{r) > for all r > 0. 

2. (p'{r) ~ ip'{s) and (p"{r) ~ <f"{s) for0<s<r<2s. 

3. |v9(i+*=)(r)| < 99'(r)/r'= for all r > and A: e N. 

Let x(r) be a dyadic cut-off function which is supported around r ^ R and satisfies 

\x^^\r)\ < R-\ (2.7) 

(These estimates are supposed to hold uniformly for r and R, but may depend on k.) Then 
we have for any d G N, 



sup 



xm)e 



di 



<t-'^'\^'{R)/R)'^''-'y\^"{R))-^l\ 



(2. 



where the constant depends only on ip, x-, d. 
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Proof. Since the case d = 1 is easy, we may assume d >2. The phase can be stationary 
when its gradient 

^ := tip' {r)e + X (2.9) 

vanishes, where we use the polar coordinates ^ = rO. Because all other parts are radially 
symmetric, we may assume that x = \x\ei, where ei = (1, 0, . . . , 0) G W^. By the monotonic- 
ity of (p', the stationary point for ^ is uniquely determined for any fixed {t,x), by ^ = cei 
where c is the unique solution (if it exists) for the equation 

|x|=V(c). (2.10) 

First we dispose of the region away from the stationary point. There exists a smooth 
partition of unity on 5'^"^ denoted by Fi{9), . . . Fn{0) G C°°{S^~'^) where N depends only 
on d, satisfying the following. 

1. J2j Fj{e) = 1 for any 9 G S'^-'^. 

2. There exists G S"^^^ for each j such that ■ > 1/2 for any 9 G suppFj. 

3. For any j > 1, ■ x < 0. 

Thus we have suppFj C {9 ■ ei < ^/3/2} for j > 1. Denote by I the integral to be 
estimated. We decompose the integral by using Fj and estimate those terms for j > 1 by 
partial integration in the direction. 

I = J2 I3 = I x(r)F, (0)e-*'^W+^«-d^ (2.11) 



Ij = J e'^{djipjYixir)Fj{9))d^, 

for any K > 0, where we denote 

P ■= -t^{r) +^-x, dj = Vj ■ V, pj = djp, 



(2.12) 



(2.13) 



and djipj ^ denotes the operator / >—>■ dj{ipj ^/), not the derivative of ipj ^. This convention 
applies to similar expressions of the form d:^p~^ appearing later, but otherwise, we use 5° 
to denote the usual partial derivatives. 
By the property of Fj and , wc have 



Pj = ■ {-tip'{r)9 + x) < -tip' {r) 12, 
for any 9 G suppFj. Expanding the X-th power, we estimate 



\a.\=K 



^''^'"''^''^' df%x{r)F,{9)) 



P 



K+l 



(2.14) 



(2.15) 
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By using the estimates < y /r% \x^^'^\ < r'^ and 1^=^^(0)1 < r-^ we obtain 

|/,| < {fy{R)R)-^R'' (2.16) 

for any j > 1 and K >Q. Choosing K = d/2, we obtain the desired estimate for those /j, 
since \(p"\ < /r\. 

Now we are left with the stationary part I\. We further decompose this region into 
parallelepiped regions by a smooth partition of unity. Their side length should be 5 := 
{t(p" (R))^^/"^ in the ei direction and A := {tip' {R) / R)^^/'^ in the other directions. Denote 
the set of such parallelepipeds by V. For each P £ V, we have a smooth cut-off function 
Gp{^). We have ^^'(0 = 1 on suppFi(6l)x(r), 6 ~ |^| ~ on any P eV, and 

\D"Gp\ < <5-"iA-"2— (2.17) 

Denote the integral corresponding to Gp by Ip. 

We put the stationary point at the center of one parallelepiped Pq. In that region, the 
integral is not very oscillatory, so we just bound it by the volume: 

\IPo\^SX^~\ (2-18) 

which yields the desired estimate for this part. 

We have the same volume bound on the other parts, but we need additional decay factors 
to sum those numerous pieces up. For that purpose we may integrate by parts either in the 
ei direction or in the other _L ei directions. 

First we consider integration in direction. Take any P E V on which we have 
l^^l ^ Xk, k e N. Notice that such a P exists only if A < /?. Then wc have a direction 
V J- ei such that v ■ ^ ^ Xk foi any ^ € P. We integrate by parts in the v direction. Then 
we obtain 

Ip = I e-'"^^^^+'^--id,ip-Yixir)Fii9)GpiO)dC, (2.19) 

for any K >0, where dy = v -V and 

= dy{-tip{r) + ^-x) = -tv C(p'{r)/r tXkip'{R)/R. (2.20) 

By expanding the K-th power, we obtain 

\Ip\ < {tXk(p'{R)Rmm{X,R))-^X'^-^S, (2.21) 

where A in min(A, i?) is coming from derivatives hitting Gp(^). Thus we get 

|/p| < k-^'dX^'-K (2.22) 

Next we consider integration in the ei direction. Here we should be a bit more careful, 
since the phase can be very small depending on (p" even if we are not in Pq, and then the 
estimate on the higher derivatives of the phase becomes less trivial. Take any P G "P on 
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which we have — c ~ 6j with j ^ 0, where c denotes the value of the stationary point. 
Such a P exists only ii 5 < R. Let pi = dip = —tip'{r)6i + We compute 



5ir = Oi, diOi = (1 - el)/r, 
dm = -t{^"Oi + ^'{1 - 9l)/r) 



2, , . (2-23) 



For higher derivatives we have 

k 

l^t+Vl < Y.^W^^^'Vr'\+ty{l - el)/r''+\ (2.24) 

Since 99', ip" > 0, we have < \dipi\/r^ . This in particular implies that dipi is of the 

same size on suppx- Then it implies that \pi\ ~ \dipi5j\. We have also \dipi\ > t\ip"{R)\. 
Integrating by parts in the ei direction we have 

Ip = I e-''^^'^+'^-^{ditp^Y{x{r)Fi{e)Gp{i))di , (2.25) 

Expanding the i^-th power, we estimate 

\Ip\ < {t^"{R)5j5)-''\''-H, (2.26) 

Thus we obtain 

\Ip\ < r''X''''S. (2.27) 
In conclusion, we have for any P G V, 

\Ip\<il + k + j)-''X''-'5, (2.28) 

for any K (^N, so we can sum them up for all k,j, deriving the desired estimate. □ 



3 Low frequency normal form 

Let u be a solution of H1.5() and v = V~^u. By using 1)2.4(1 . we can derive the equation for 
v: ^ 

y{t)=e-''^vo- f e-'^'-'^"V-hF{u{s))ds. (3.1) 
Jo 

The presence of in the nonlinearity is our only "enemy" for d > 4; without it, the 
scattering would follow in the same way as for NLS with quadratic and cubic power nonlin- 
earities. Quadratic nonlinearities have the critical decay rate for scattering by the Strichartz 
estimate when d = 4, while for d > 4 their decay is more than sufficient. Indeed we can 
derive the scattering result for d > 4 and certain restricted values of a by using only the 
Strichartz estimate with the U gain and the Sobolev embedding. That argument, however, 
cannot work in the critical case d = 4. 
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We introduce the following nonlinear change of variable to kill the singular nonlinearity. 
It gives a better range of a even for d > 4. Let 

w:= N{u) := u + P\u\'^/2, (3.2) 

where P is the Fourier multiplier F^^x^i X S Cq°(M'^) satisfying x(^) = 1 for \x\ < 1 and 
x{x) = for |x| > 2. Then w satisfies the following equation: 

iw = -Aw + 2Rew + Gi+iG2, 

Gi = (3 - P)ul + Qui + PA\u\'^/2 + l-upm, (3.3) 
G2 = 2Q{uiU2) + VP • {U2VU1 - U1VU2) + Qdtipna), 

where u = ui + iu2 and Q = Id — P. Thus 

z ■= V-^w = V-^N{u) = U-^iui + P\uf/2) +iu2 = v + U-^P\Vvf/2 =: Mv 

satisfies ^ 

z = e-^^'z{0)- [ e-'^^'-'\iGi-U-^G2)ds. (3.4) 

Notice that the singularity of U^^ in front of G2 is now canceled by either Q or V. There is 
no loss of regularity thanks to P. Moreover, G does not contain bilinear interaction of very 
low frequency of U2, whereas ui behaves better at the Fourier due to the relation with z. 



4 Proof of the main theorem 

We introduce Besov spaces with different regularity for low and high frequency: 

IIV'lls^^,*. := WP^Wb^ + \\Qv>\\b^^, H'^^' := (4.1) 

where B^ := 3*^2 denotes the homogeneous Besov space (cf. T). We will use these 
spaces under the restriction q < d/a, embedding them into S'. We have Bg'^ = U^~°'i3^2i 
jjOfi _ j^b ^j^g embedding Bg'^ C Bq for a < a' and b > b'. These spaces have the 
following multiplicative property, which is enough for us to estimate the nonlinearity. 

Lemma 4.1 Let pi,p2,P3 G [2, 00], 81,82,83 G R, ti,t2,t3 G R and ba := 1/pa for a = 
1, 2, 3. As8ume that for a = 1,2, 3, 

max(0, 8a, 81 + 82 + S3) < d{bi + 62 + ^3 - 1) < + ^2 + ^3 > ta, 8a < dba- (4.2) 

Then we have 

I fghdx 



< 



rS1.*1 



^P2 



(4.3) 
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Proof. We apply the Littlewood-Paley decomposition to each of the functions. 

f9h= Yl Yl H ^-J^'' / f^9khi + fkgihj + fWjhk, (4.4) 

l=-ook=l-2j=-oo •' 

where < Cj^^^i < 1 and each fj is supported in Fourier space in {2^^^ < |^| < 2^~^^}. By 
symmetry, it suffices to estimate the first summand. Let /? = 1 — 62 — ^35 then we have 
< P < bi. By the Sobolev and the Holder inequalities and the embedding Bp C L^, we 
have 



< 



/jllLV/slbfellLPz ||/i;||if3 



j<k 



(4.5) 



for any s. Taking s = si for j < and s = ti for j > (when k > 0), the norm on / is 
bounded by 

f^ki-si+dih-m {k<(}), 

jl + 2M-ti+rf(6i-/3)) (fc>0), ^ 
since — si + d{bi — /3) > 0. The summation over k ^ I e Z converges, because for A; < 0, 

-si + d{bi - /?) - S2 - S3 = -{si + S2 + S3) + d{bi + 62 + &3 - 1) > 0. (4.7) 

and for A; > 0, 

-t2 - is < 0, -h -t2-h + d{bi + 62 + - 1) < 0, (4.8) 

and the summability is provided by the £^ norm in the Bcsov norms. □ 

For the proof of our main theorem, we will use the following special cases of the above 
lemma. Let d > 4, s = d/2 - 1, b = 1/p = 1 - 1/p' = 1/2 - 1/d and 1/q = 1/2 - l/(2(i). 
Then we have 

(1) H"-^''-^ X B^-^^'-^ C B^-;^^''-^-^ , 

(2) H"^' X B^-^^'-^ C H"-^''-^, 

(3) 5^-^/2,. ^ 50-6/2,5 ^ ^a,.^ ^4 9) 

(4) H"'' X F'^'" C Bpf, 

(5) X B^-''/^'' C i/'^-i.«-i/2^ 
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for < j,k,j + k < 1, provided that \a\ — b < d/2 — 2. All conditions in the lemma are 
satisfied because 

(1) max(cr - ^h,\a\ - b) < d/2 - 2 < s > j + k - s, \a\-b < d/2 - 1, 

(2) max((T, 6-cr)<d/2-l<s>l-s, a < d/2, 

(3) max(cr - 6/2,-0-) < (i/2 - 1 < s, a - 6/2 < d/2 - 1/2, (4.10) 

(4) \a\<d/2-l<s, 

(5) max(CT - 6/2 + l,cj,l - a) < d/2 - 1/2 <s + l/2> 1/2 - s. 

As the function space for v = V~^u and z = V~^N{u), we define Xq and Xi by 

Xj = n (411) 

We have U : Xq Xi C Xq bounded. First we consider the norm relation between 



V and z. By (Ilini)-(2)(4) and the embeddings B'';' C i?'^'" C Bp^'''' ^ and 



jja~b,s~i ^ Bp'^'"''^, we have 



|f^-'^'(/9)l|B < (4.12) 



for i3 = i?""'* and i? = Bp'^ ''^ . Hence for any small z G H'^'^ , the unique inverse v = M~^z 
is obtained in the same space by applying the Banach fixed point theorem to the equation 
V = z — PlVvl"^ /2. Thus M is bi-Lipschitzian for small data in H'^'^. In particular, we 
have the following estimate for the free part zq := e~*^*z(0) 

ll^olUo < MO)\\h^,s < \\vm\H^,s < 6. (4.13) 

Moreover, 1)4.12(1 implies that ||z||xo ~ Ibllxo ^ long as they remain small. In addition, we 
have 

llM(e-^-^V) - e-^^Vlk-.^ ^ (|t| ^ oo) (4.14) 

for any small ip £ H'^'^ . This is proved as follows. First ip £ C^, then the decay estimate 
and (|4.9|l -(5l implies that 

||t/-ip|ye-^^Vl'llH^- < ||e-^^Vlk-l|e-'''VlL-'>A. < t-'^'- (4.15) 

Then (|4.14() for general if follows from the denseness of C H'^'^ and the Lipschitz 
continuity of M. The bi-Lipschitz property implies also that 

We-'^'if - M-^(e-*^V)lk-.» ^ (|t| ^ oo) (4.16) 

for any small if G H°''^. 

Next we estimate the quadratic terms in G{u). By (|4.9() -fl'). we have 



Ulu\\r2n<T+b,s < ||^l||L°°//<^-''.s|l^llr2D<^-i'.s ^ ll^lllXi ll^llxo) 



p 



||P(nVn)||^2ga+6,s < \\u\\^2^^-b,s\\Vu\\Lo.H.-b,s-i < \\ufx^, 



(4.17) 
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These are used to estimate the terms uf, uiU2, PA|np and U2VU1 — uiVu2, together with 
the boundedness of U'^Q, U'^VP and PV on every Besov space. can be treated in 
the same way as uiu since 

Qui = Q{U2 ■ Q^2U2) + Q{Q-2U2 ' P-2U2) (4.18) 

and Q-2 : -'^o — ^ bounded, where P-2 := J^^^xi^O^ ^^'^ Q-2 ■= Id — P-2, defined in 
the same way as P and Q (see below p.2() ). 

For the cubic terms, we use (UH)-(3)(4) together with the interpolation inequality Xq C 
L^Bq yielding 

II^^IIl2b-;- < WnWi^wAW^Wh^H"-- ^ lkllxolkll^4s--i'/2,» < (4.19) 

Gathering these estimates, we obtain 

Ik - zoWxo < \\V-^iG\\^^^a+b,s < ||n||x„(||n||xo + IhilUi + lhllx„) ^ lkllx„> (4.20) 
if the right hand side is small. By continuity in time and H4.1l-{jl . we conclude that 

\\z\\xo + ll^olUo + \\V-'^iG\\^,^.+b,s < \\v{0)\\h-,s. (4.21) 

p' 

In particular, the norm of G on the interval (T, oo) vanishes as T — > oo, which implies via 
the Strichartz estimate that e^^^z{t) — > ^v± strongly in H"''^ . Then 1)4. 16() implies that 
e'^^v{t) v± in H''^' . 

We can construct the wave operators applying the same estimates to the Cauchy problem 
from t = ibcxD: 

z = e-'"'v± - f e-'^^'-'W-\G{VM-^z{s))ds. (4.22) 

J itoo 

The bi-Lipschitz property of the maps w(0) i— > v± is similarly proved by estimating the 
difference of two solutions in the same space. These arguments, essentially by the Banach 
fixed point theorem, are quite standard in nonlinear scattering theory, and so we omit the 
details. 
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